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Abstract 

The event-by-event analysis of heavy ions collisions is becoming possible 
with advent of large acceptance detectors: it can provide dynamical infor- 
mation which cannot be obtained from inclusive spectra. We identify some 
observables which can be related to standard thermodynamical theory of fluc- 
tuations and therefore may provide information about properties of hadronic 
matter at the freeze-out stage. Among those are fluctuations of "apparent 
temperature" (the p t slope), as well as the population of different bins of the 
pion momentum distribution. 
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l.In the field of heavy ion collisions the event-by-event search for "un- 
usual" events (deviating in some way from the average behavior) is attracting 
increasing attention. Some future experiments at RHIC (especially STAR 
and PHOBOS) consider it to be among the central issues of their physics 
program. 

In general, deviations from the average behavior can be of dynamical or 
of statistical nature. A well known example of the former kind is strong 
dependence of all observables on the value of the impact parameter b|. Fur- 
thermore, direction of the impact parameter b in the angle is important for 
studies of the specific collective effect known as "directed flow". Non-trivial 
deviations of dynamical nature due to QCD phase transition were for ex- 
ample predicted in the AGS/SPS energy domain, due to formation of QGP 
bubbles [Q or of the disoriented chiral condensate (DCC) 0. In order to see 
whether those are present or absent in the data, one should have a proper 
benchmark, the statistical fluctuations we discuss below. 

Fluctuations can also be sensitive to non-trivial details of the dynamics. 
An interesting example was considered by Gazdzicki and Mrowczynski ||: 
event-by-event fluctuations of the < p t > e vent turns out to be sensitive to 
models describing nuclear collisions. In particular, models considering just a 
superposition of independent pp events (like e.g. the original VENUS event 
generator without re-scattering) predict much larger fluctuations compared 
to models with significant re-scattering of secondaries (e.g. RQMD). This pa- 
per has triggered experimental studies, and the preliminary results coming 
from the NA49 experiment H] have found that the fluctuations (in the appro- 
priate variable) follow perfect Gaussian | with small width, clearly excluding 
any model based on superposition of independent pp collisions. 

This impressive example have initiated the present work, which also pro- 
poses several observables for which analysis of fluctuations is potentially use- 

2 Because of its paramount significance most experiments are able to fix its value 
already at the trigger level, by measuring charge multiplicity, transverse or "spectator" 
energy. 

3 No non-statistical "tails" deviating from it was found, although Gaussian was followed 
for several orders of magnitude. 
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ful for uncovering the details of the system's evolution. However, our ap- 
proach is more general, it is not based on any particular model but rather on 
well known theory of thermodynamical fluctuations. Although the predicted 
fluctuations are probably not so spectacular as some dynamical models sug- 
gest, they are guaranteed to be there. 

2. In general, any statistical fluctuations can be derived from the famous 
Boltzmann expression relating entropy S and probability P, rewritten by 
Einstein in the following form 



P ~ exp(S) (1) 



Before we turn to specifics, few general comments about it are appropri- 
ate. Applying general thermodynamical theory of fluctuations we rely on 
statistical independence of fluctuations in different volume elements. The 
fluctuations of temperature or particle composition should be independent 
on many dynamical details: it is irrelevant whether all elements of the excited 
system have their freeze-out at the same time or not, whether they freeze-out 
at rest in the same coordinate frame or are all moving with different veloci- 
ties^] The only thing which matters is how the global entropy of the system 
S depends on the particular observable under consideration. 

In general, the fluctuations predicted by (1) neither should be small, nor 
the discussed system should have macroscopic number of degrees of freedom. 
The only requirement is that the system is equilibrated, in the sense that 
it equally populates all its available phase space. Applications of theory of 
thermodynamical fluctuations to multi-hadron production reaction have a 
long history, going for pp and pp reactions back to 70's. In || probability 
of many exclusive channels for low energy pp annihilation were calculated: 
simple textbook formulae for the entropy have predicted quite accurately 

4 Still it is useful to imagine that the system is cut into many small (but still macro- 
scopic) elements, and if each is taken when it has a given T, \x and put to rest and combined 
into a common system, the whole system will occupy some total volume V. 
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their probabilities, starting with channels with only 4 secondaries^ In || 
the production of KK pairs in a pion gas was considered along these lines 
further. Summing thermodynamical formulae over states with only proper 
quantum numbers (no strangeness), a general form of a general dependence 
of the K/tt ratio on the system size (the pion multiplicity N n ) was derived. It 
changes from 0(N n ) value for small enough systems to the A^-independent 
ratio for large ones, and describes the data perfectly. Very recently similar 
approach was used in ref . |J , where excellent thermodynamical description of 
the yields of impressive list of secondaries in e + e~,pp and nuclear collisions 
has been demonstrated^. 

3. Before discussing statistical fluctuations in details, let us make a general 
comment about their size. For heavy ion reactions in question the multiplicity 
is iV = O(10 2 — 10 4 ). Therefore, the total entropy produced is large and the 
relative statistical fluctuations are small. Obviously they scale as 1/ and 
are expected to be at few percent level. Sufficient experimental statistics to 
measure such fluctuations is not a problem. 

In a standard way, the mean square deviations of various quantities can 
be obtained from the expansion of the entropy in powers of the deviation up 
to second order: thus proper succeptibilities appear. So, if measured, the 
coefficients of 1/yN can provide useful information about thermodynamics 
of hadronic matter. 

Our first example is the fluctuations in temperature, which was recently 



5 Recall that the phase space for 4 particles is already a 8-dimensional integral. Us- 
ing thermodynamical expressions instead of exact the phase space means saddle point 
approximation: it has about the same accuracy as Sterling formula for the factorial 

Nlog(N/e), which is excellent even for not-so-large N. 

6 We still do not understand why e + e~ ,pp look so equilibrated (see e.g. |l(J for recent 
discussion) . Heavy ion collisions is a different matter: they produce larger and longer-lived 
systems, and also display other features of a macroscopic behavior such as the transverse 
flow. 
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discussed in ref.0[]. Standard thermodynamics (see e.g. ||) tells us that 



exp[-^(AT/Tf] (2) 



where AT is the deviation of temperature from its mean value and C V (T) is 
the heat capacity of a hadronic matter. (The heat capacity which enters here 
is extensive quantity C v = T^\n,v, proportional to the total volume V or 
the number of particles N in the system, so the relative fluctuations depend 
on N as expected.) So in general, by measuring fluctuations in T we may 
learn about the heat capacity. 



Let us now look at this general idea in more details, discussing which 
temperature we actually measure and what behavior of the fluctuations we 
may expect. If the temperature is deduced from slopes of the observed 
m t = ^F^Tm? spectra, it corresponds to the so called thermal freeze-out 
conditions^. It happens when the expansion rate of the system becomes com- 
parable to re-scattering rate. Note however, that for different secondaries the 
re-scattering rate is rather different: say for SPS conditions it is larger for 
nucleons than for pions, and is even smaller for kaons. Therefore one should 
see different temperatures, growing from iV to tt to K. This is indeed what 
happens if one study realistic cascades (e.g.RQMD). 

Recently detailed study of thermal freeze-out have been also made in 
hydro-based framework [pl|. In contrast to naive freeze-out conditions used 
in most earlier hydro works (freeze-out surfaces approximated by isoterms 
T = Tj) proper kinetic condition |13|] (see also [|14], [15[]) was used, separately 



for 7T, K, N. It was found that for different secondaries and collisions those 
surfaces have quite different shapes. In particular, the system cools deeper at 
its center, and remains hotter at edges: so the relevant final Tj is by no means 

I apologize that I did not know about it while preparing the preprint form of this 
work: let me thank S. Mrowczynski who informed me. 

8 Not to be confused with those for chemical freeze-out ones at which particle composi- 
tion is determined. The former is related to elastic scattering rates, the latter to inelastic 
species-changing ones. 
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a universal constant. Furthermore, they especially depend on the absolute 
system size: those which are produced with AA collisions with larger A cools 
to lower Tf. For example, for medium-size nuclei (Si-S) these results suggest 
Tf = 140 — 150MeV, while for central heavy ion collisions (AuAu at AGS 
and PbPb at SPS) the hadronic re-scattering goes on for long time (about 
10-15 fm/c) and typical Tf becomes as low as 100-120 MeV. Additional test 
of these numbers was provided by a detailed study of the radial collective 
flow IIS] and analysis of HBT radii pO] coming from NA49 experiment. 



Let us now return to fluctuations. Experimentally one should better look 
at some simple observable, e.g. deviations from the mean p t or a fitted 
slope of the mt distribution T in a event. As it is well known, T is not 
temperature but a combination of freeze-out local temperature Tf and trans- 
verse flow velocity v t of matter elements, averaged over the whole system. 
It is complicated in practical calculations, but we think the actual situation 
is significantly simplified because the fluctuations of v t are likely to be less 
important then those ofTf. This conjecture is motivated as follows: (i) Un- 
like Tf (which depends on random statistical separation of total energy into 
transverse momenta, longitudinal momenta and particle masses at the freeze- 
out moment), v t is basically an accumulated acceleration during the whole 
evolution. Random fluctuations in pressure should be significantly washed 
out after the time integration, (ii) Furthermore, at AGS/SPS energies most 
of the time the expanding system spends near the QCD phase transition, in 
the so called mixed phase. As a result, most of the time available the pressure 
is nearly constant and cannot significantly fluctuate, (iii) The relative role 
of Tf and v t in the observed slope significantly depends on the particle mass: 
Vt is much less important for pions than for nucleons or deutrons. One may 
compare fluctuations of ir, N slopes (the main secondaries) to test our con- 
jecture, (iv) Finally, for pions corrections for transverse motion is basically 
a "blue shift" factor, which cancels in ratios. 



This conjecture implies that 



f 1 f f ' C v {T f ) { ] 



The key observation at this point is that because C V (T) has strong In- 
dependence, it can be used as a "thermometer" to test matter properties 
at thermal freeze-out. 



Particular form of the T-dependence of the heat capacity in the vicinity 
of the QCD phase transition (T ~ T c m 150MeV) depends on the order of 
the phase transitionP]. Large peak at T = T c has been observed in lattice 
numerical data 0], but its exact forrri^] and extrapolation to real QCD is 
not yet clear. In the hadronic phase T < T c the thermodynamics is usually 
described as a "resonance gas" |lTJ. As shown in fL8f , the corresponding EOS 



agrees well with what is used in cascades like RQMD, and also explains the 
observed flow. For practical purposes its T-dependence can be conveniently 
parameterized as e,p ~ j ll + 1 / c where parameter c 2 (square of the sound 
velocity) is (nearly) T-independent. It still depends on matter composition 
though (baryon charge/entropy ratio), changing from c 2 ~ .14 for baryon- 
rich matter produced at AGS to c 2 ~ .2 for baryon-poor one at SPS. Such 
EOS leads to C v ~ T 1 ^ 2 ~ T 7 (T 5 ) at AGS (SPS), respectively: so even in 
the hadronic phase the T-dependence of C v is strong enough to make a good 
thermometer. 



Let us now discuss how it is supposed to work. At qualitative level, a 
debate continues on whether the observed pions are emitted directly from the 
QGP clusters/mixed phase (see e.g. |TJ, |T7|) at T ps T c , or (as we advocated 
above) they cool further into a resonance gas phase (especially for heavy 
ions, Pb or Au). In the former case one should find C V (T) at its peak: the 
temperature fluctuations should then be strongly suppressed. In the latter 
scenario, C V (T) is much smaller, and we end up with a (counter- intuitive!) 
prediction that fluctuations for heavy ions should be relatively^ stronger 
compared to those for medium ions. 

Using the EOS of the resonance gas and numbers for freeze-out tempera- 

9 Dclta function for the 1-st order, divergent peak for the second order and a finite peak 
for a rapid cross-over. 

10 It depends on parameters like quark masses used and other methodical details like 
Wilson versus Kogut-Susskind fermions used. 

11 Recall that we are discussing a coefficient of 1/V~N. 
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tures Tf for medium and heavy ion collisions mentioned above ( ||18|| ) we can 
make more quantitative predictions. 



[N < (ATf/Tf) 2 >] hea vy C v (Tf edAum ) Tf edium /c2 

[N < (ATf/Tf) 2 >Y nedium ^ C (Tf eaVy ) ^ rpheavu 



(4) 



Although the values of Tf mentioned above are not that different, this ratio 
should change by about factor 2, from medium to heavy ions. Such drastic 
change should be easily observable. 

4. Now we proceed to fluctuations of a more detailed observable, the 
occupation of particular momentum bins in the histograms. In purely statis- 
tical system the fluctuations in the particle number are generally described 
by 



< AN 2 >=T^\ T ,, (5) 



dfx 



For example, for classical ideal gas we got trivial Poisson statistics. Devia- 
tions from it can be induced without interaction, just by quantum statistics. 
For the ideal Bose-Einstein (Fermi-Dirac) gas the expression above leads to 
enhanced (suppressed) fluctuations: 



= (1± < n k >) (6) 

< n k > 



So, by measuring deviations from Poisson statistics (1 in the r.h.s) one can 
learn what is the mean quantum degeneracy factor < n& >. 

Bose-Einstein (Fermi-Dirac) correlations in two-body distributions of pi- 
ons (nucleons) originating from heavy ion collisions have been studied for a 
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long time and is the basis of pion (nucleon) interferometry providing infor- 
mation about space-time picture of the source. Correlation of 3 and more 
particles should be stronger but present in so small corner of the phase space 
that those have not been used so far. 

Fluctuations in the distributions on event-by-event basis we discuss is 
just another way to approach the same "induced radiation" phenomenon. 
The more bosons in a given momentum cell is produced in a given event, 
the larger is the probability to find another one in the same cell. Like for 
pion pairs, the interference happens late and the role of interaction should 
be small. 

One of the problems with practical implementation of this simple idea 
is that the total multiplicity of produced pions is known to be a subject of 
large non-statistical fluctuations, in the first place due to different impact 
parameters in different events. There may be also other dynamical effects 0- 
Therefore one should consider deviations from the normalized distributions, 
or relative fluctuations in different bins. 

Pions are the only secondaries which are light enough to have noticeable 
quantum degeneracy n k at freeze-out stage. For chemically equilibrated pion 
gas (chemical potential /x w = 0) at zero momentum Bose-enhancement factor 
1 + nj, is increasing statistical fluctuations (relative to random statistics) by 
a factor of about 1.5. Naturally, the effect disappears for large momenta. 

Furthermore, it was argued that the pion gas is not chemically equili- 
brated at thermal freeze-out, but has a non-zero chemical potential^ so 

12 The multiplicity distribution in pp collisions are known to be wide and very non- 
statistical, and there were even speculations that they display fractal (intermittent) dis- 
tributions. However, this is not so for nuclear collisions: the only noticeable correlations 
observed so far in this case are due to Bose-Einstein correlations we discuss, see e.g[[l6| . 

13 The issue has a somewhat controversial history. It first surfaced because of experi- 
mentally observed "low p t enhancement" relative to equilibrium = 0) distributions. 
For discussion of pion "over-population" of the phase space see e.g. |H| and references 
therein: even a "pion laser" conditions ( n&=o — > oo) were discussed. However, in fact 
most of this enhancement is due to resonance decays. 
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that the effect should be even stronger. The basic reason for non-zero (n n was 
originally pointed out by G.Baym: adiabatic expansion of the re-scattering 
pion gas creates it, even if one starts from chemical equilibrium (fi n = 0). 
Estimates made in |21], |18[ has resulted in small effect for medium ions, but 
for heavy ones (PbPb collisions at SPS) it should reach for mid-rapidities 



fi n ~ 60MeV. Moreover, in |L8| it is shown that comparison of NA44 data 
for PbPb and SS shows this effect, with fi n ~ 60MeV. If so, Bose enhance- 
ment of fluctuations at small p t should be larger, giving the 1 + factor of 
2.0-2.3. Again this prediction looks counter-intuitive: (relative) fluctuations 
are enhanced for larger (not smaller) systems! Finally, let us point out direct 
relation between this effect and one of the possible manifestation of DCC 
As discussed in ||, the remnants of the disoriented vacuum can be vied as 
the unusually large population of the \cw-p t bins in the pion p t distribution. 
If this indeed happens, expression ([)]) tells us that one should also observe 
enhanced fluctuations in the corresponding bins induced by Bose-Einstein 
interference effects. 



So far we have discussed an idealized picture of a gas of secondary pions: 
now we are going to face real situation. First of all, discussion above implied 
that there exist some common reference frame in which it is at rest, so that 
the "small momentum" bin discussed above in unambiguous. In reality the 
pion system consists of volume elements moving with different collective ve- 
locities, and this restrict/dilute the Bose enhancement effects. The density 
leading to enhancement is actually density in the phase space: even for pi- 
ons with identical momenta this effect is absent if their emission points are 
too far apart in space or too much separated in time. This issue was much 
debated for two-pion interferometry, with realization that the the so called 
HBT radii R x ,y,z are not the size of the system, but a size of a region in which 
the velocity change (due to its gradient) can be compensated by thermal ve- 
locities. (In other words, those are related to curvatures of the freeze-out 
surface.) Bose enhancement of fluctuations is as strong as discussed above, 
provided the bin used has the size such that \Api\Ri < 1. If for statistics 
reason the bin is larger, the effect is reduced accordingly. 

The second problem is again well understood for two-pion interferometry. 
Substantial fraction f LL of secondary pions come from long-lived (LL) reso- 
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nances like u and weak Kq decays: those do not interfere with others (unless 
the resolution or bin size is as small as the inverse lifetime). Two-pion inter- 
ferometry and the Bose enhancement term we speak about is reduced by it 
by the same factor (1 — /ll) 2 - 

In summary, we have pointed out that statistical fluctuations of the col- 
lective observables are related to thermodynamically interesting quantities. 
In particular, the r.m.s. deviation in the apparent temperature of the events 
is related to the heat capacity at freeze-out, and fluctuations in the p t spec- 
trum at small p t are sensitive to freeze-out temperature and pion chemical 
potential. 
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